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We report on a measurement of CP-violating asymmetries (Acp) in the Cabibbo-suppressed 
D" — ^ ■K'^n~ and — >■ K'^ K~ decays reconstructed in a data sample corresponding to 5.9 fb~^ 
of integrated luminosity collected by the upgraded Collider Detector at Fermilab. We use the 
strong decay D*"*" — >• _D°7r+ to identify the flavor of the charmed meson at production and exploit 
CP-conserving strong cc pair-production in pp collisions. High-statistics samples of Cabibbo-favored 
P)° — ^ K~i^^ decays with and without a D*^ tag are used to correct for instrumental effects and sig- 
nificantly reduce systematic uncertainties. We measure Acp{D° — >■ 7r"'"7r~) — (+0.22 ±0.24 (stat) ± 
0.11 (syst))% and Acp{D° K+ K") = (-0.24 ± 0.22 (stat) ± 0.09 (syst))%, in agreement with 
CP conservation. These are the most precise determinations from a single experiment to date. Un- 
der the assumption of negligible direct CP violation in D^' — >■ tt+tt" and D^' — >■ decays, the 
results provide an upper limit to the CP-violating asymmetry in mixing, |App(D'')| < 0.13% at 
the 90% confidence level. 



I. INTRODUCTION 

The rich phenomenology of neutral flavored mesons 
provides many experimentally accessible observables sen- 
sitive to virtual contributions of non-standard model 
(SM) particles or couplings. Presence of non-SM physics 
may alter the expected decay or flavor- mixing rates, 
or introduce additional sources of CP violation besides 
the Cabibbo-Kobayashi-Maskawa (CKM) phase. The 
physics of neutral kaons and bottom mesons has been 
mostly explored in dedicated experiments using kaon 
beams and e+e^ collisions The physics of bottom- 
strange mesons is currently being studied in detail in 
hadron collisions In spite of the success of several 
dedicated experiments in the 1980's and 1990's, experi- 
mental sensitivities to parameters related to mixing and 
CP violation in the charm sector were still orders of mag- 
nitude from most SM and non-SM expectations 0. Im- 
provements from early measurements at dedicated e^e~ 
colliders at the T(4S') resonance (S-factories) and the 
Tevatron were still insufficient for discriminating among 
SM and non-SM scenarios 0-01 ■ Since charm transitions 
are described by physics of the first two quark gener- 
ations, CP-violating effects are expected to be smaller 
than 0(10"^). Thus, relevant measurements require 



large event samples and careful control of systematic un- 
certainties to reach the needed sensitivity. Also, CP- 
violating effects for charm have significantly more un- 
certain predictions compared to the bottom and strange 
sectors because of the intermediate value of the charm 
quark mass (too light for factorization of hadronic am- 
plitudes and too heavy for applying chiral symmetry). 
All these things taken together have made the advances 
in the charm sector slower. 

Studies of CP violation in charm decays provide a 
unique probe for new physics. The neutral D system 
is the only one where up-sector quarks are involved in 
the initial state. Thus it probes scenarios where up- type 
quarks play a special role, such as supersymmetric mod- 
els where the down quark and the squark mass matrices 
are aligned Q and, more generally, models in which 
CKM mixing is generated in the up-quark sector. The in- 
terest in charm dynamics has increased recently with the 
observation of charm oscillations [lol - [l^ . The current 
measurements [3 indicate 0(10-^) magnitudes for the 
parameters governing their phenomenology. Such values 
are on the upper end of most theory predictions p^ . 
Charm oscillations could be enhanced by a broad class 
of non-SM physics processes Any generic non-SM 

contribution to the mixing would naturally carry ad- 
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ditional CP-violating phases, which could enhance the 
observed CP-violating asymmetries relative to SM pre- 
dictions. Time integrated CP-violating asymmetries of 
singly-Cabibbo-suppressed decays into CP eigenstates 
such as I?° —J' TT+TT" and K^K^ are powerful 

probes of non-SM physics contributions in the "mixing" 
transition amplitudes. They also probe the magnitude 
of "penguin" contributions, which are negligible in the 
SM, but could be greatly enhanced by the exchange of 
additional non-SM particles. Both phenomena would, 
in general, increase the size of the observed CP violation 
with respect to the SM expectation. Any significant CP- 
violating asymmetry above the 10 ~^ level expected in the 
CKM hierarchy would indicate non-SM physics. The cur- 
rent experimental status is summarized in Table HI No 
CP violation has been found within the precision of about 
0.5% attained by the Belle and BMAR experiments. The 
previous CDF result dates from 2005 and was obtained 
using data from only 123 pb^^ of integrated luminos- 
ity. Currently, CDF has the world's largest samples of 
exclusive charm meson decays in charged final states, 
with competitive signal purities, owing to the good per- 
formance of the trigger for displaced tracks. With the 
current sample CDF can achieve a sensitivity that allows 
probing more extensive portions of the space of non-SM 
physics parameters. 

We present measurements of time-integrated CP- 
violating asymmetries in the Cabibbo-suppressed — > 
K^K~ decays (collectively referred to 

I 



as — >■ h~^h~ in this article) using 1.96 TeV proton- 
antiproton collision data collected by the upgraded Col- 
lider Detector at Fermilab (CDF II) and corresponding 
to 5.9 fb"-'^ of integrated luminosity. Because the final 
states are common to charm and anti-charm meson de- 
cays, the time-dependent asymmetry between decays of 
states identified as and at the time of production 
{t = 0) defined as 



Acpih+h-,t) = 



N{D° h+h-;t)-N{D° ^ h+h-;t) 
N{D" ^ h+h-;t) + N(DO h+h--t)' 



receives contributions from any difference in decay widths 
between D° and D° mesons in the chosen final state (di- 
rect CP violation) , any difference in mixing probabilities 
between and Z?" mesons, and the interference between 
direct decays and decays preceded by flavor oscillations 
(both indirect CP violation). Due to the slow mixing 
rate of charm mesons, the time-dependent asymmetry is 
approximated at first order as the sum of two terms. 



Acp{h+h-;t) « A'^jp{h+h 



T 



A'SUh+h-), (1) 



where t/r is the proper decay time in units of D'^ lifetime 
(r « 0.4 ps), and the asymmetries are related to the 
decay amplitude A and the usual parameters used to 
describe fiavored-meson mixing x, y, p, and g by 
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where rjcp = +1 is the CP-parity of the decay final state 
and if is the CP-violating phase. The time-integrated 
asymmetry is then the time integral of Eq. ([T|) over the 
observed distribution of proper decay time {D{t)), 

Acp{h+h~) = Afp{h+h-) + A'^${h+h-) / - D{t)dt 

Jo T 

= Afp{h+h-) + ^ A'SUh+h-). (4) 

T 

The first term arises from direct and the second one from 
indirect CP violation. Since the value of (t) depends on 
-D(t), different values of time-integrated asymmetry could 
be observed in different experiments, depending on the 
detector acceptances as a function of decay time. Thus, 
each experiment may provide different sensitivity to A'^jp 
and A^^p. Since the data used in this analysis were col- 
lected with an online event selection (trigger) that im- 
poses requirements on the displacement of the P'*'-meson 



TABLE I. Summary of recent experimental measurements of 
CP-violating asymmetries. The first quoted uncertainty is 
statistical, the second uncertainty is systematic. 

Experiment Acp{tt+tt-) (%) Acp{K+K-) (%) 

B^AR 2008 -0.24 ± 0.52 ± 0.22 4-0.00 ± 0.34 ± 0.13 

Belle 2008 [16] -0.43 ± 0.52 ± 0.12 -0.43 ± 0.30 ± 0.11 
CDF 2005 [17] +1.0 ±1.3 ±0.6 +2.0 ± 1.2 ± 0.6 



decay point from the production point, our sample is 
enriched in higher-valued decay time candidates with re- 
spect to experiments at the B-factories. This makes the 
present measurement more sensitive to mixing-induced 
CP violation. In addition, combination of our results 
with those from Belle and BA3AR provides some discrim- 
ination between the two contributions to the asymmetry. 
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II. OVERVIEW 



In the present work we measure the CP-violating 
asymmetry in decays of and D mesons into tt+tt" 
and K'^K~ final states. Because the final states are 
charge-symmetric, to know whether they originate fi:om 
a D° or a decay, we need the neutral charm candidate 
to be produced in the decay of an identified or D*^ 
meson. Flavor conservation in the strong-interaction de- 
cay of the D*^ meson allows identification of the initial 
charm flavor through the sign of the charge of the tt me- 
son: D*+ D° 77+ and D*- W tt' . We refer to D 
mesons coming from identified Z?*+ decays as the tagged 
sample and to the tagging pion as the soft pion, tt^. 

In the data collected by CDF between February 2002 
and January 2010, corresponding to an integrated lumi- 
nosity of about 5.9 fb~^, we reconstruct approximately 
215 000 i:)*~tagged D" tt+tt" decays and 476 000 
D*-tagged D° — > K^K^ decays. To measure the asym- 
metry, we determine the number of detected decays of 
opposite flavor and use the fact that primary charm and 
anti-charm mesons are produced in equal numbers by the 
CP-conserving strong interaction. The observed asym- 
metry is the combination of the contributions from CP 
violation and from charge asymmetries in the detection 
efficiency between positive and negative soft pious from 
the D*^ decay. To correct for such instrumental asym- 
metries, expected to be of the order of a few 10~^, we 
use two additional event samples: 5 million tagged, and 
29 million untagged Cabibbo-favored D° — >■ K~tt^ de- 
cays. We achieve cancellation of instrumental asymme- 
tries with high accuracy and measure the CP-violating 
asymmetries of — )■ tt+tt" and — >■ K^K~ with a 
systematic uncertainty of about 10^'^. 

The paper is structured as follows. In Sec. lIIII we briefly 
describe the components of the CDF detector relevant for 
this analysis. In Sec. lIVI we summarize how the CDF trig- 
ger system was used to collect the event sample. We de- 
scribe the strategy of the analysis and how we correct for 
detector-induced asymmetries in Sec. [V] The event se- 
lection and the kinematic requirements applied to isolate 
the event samples are presented in Sec. |Vll the reweight- 
ing of kinematic distributions is discussed in Sec. IVIII 
The determination of observed asymmetries from data is 
described in Sec. I Villi In Sec. |IX] we discuss possible 
sources of systematic uncertainties and finally, in Sec. 
1x1 we present the results and compare with measure- 
ments performed by other experiments. We also show 
that by combining the present measurement with results 
from other experiments, we can partially disentangle the 
contribution of direct and indirect CP violation. A brief 
summary is presented in Sec. lXIl A mathematical deriva- 
tion of the method employed to correct for instrumental 
asymmetries is discussed in Appendix |^ and its valida- 
tion on simulated samples is summarized in Appendix 

m 



III. THE CDF II DETECTOR 

The CDF II detector has a cylindrical geometry with 
forward-backward symmetry and a tracking system in a 
1.4 T magnetic field, coaxial with the beam. The track- 
ing system is surrounded by calorimeters [isj and muon- 
detection chambers jl^ . A cylindrical coordinate system, 
(r, (/), z), is used with origin at the geometric center of the 
detector, where r is the perpendicular distance from the 
beam, (j) is the azimuthal angle, and the z vector is in 
the direction of the proton beam. The polar angle 9 with 
respect to the proton beam defines the pseudorapidity 
7y = -lntan(6l/2). 

The CDF II detector tracking system determines the 
trajectories of charged particles (tracks) and consists of 
an open cell argon-ethane gas drift chamber called the 
central outer tracker (COT) |20i] and a silicon vertex mi- 
crostrip detector (SVX II) ^ifT^The COT active volume 
covers \z\ < 155 cm from a radius of 40 to 140 cm and 
consists of 96 sense wire layers grouped into eight al- 
ternating axial and 2° stereo superlayers. To improve 
the resolution on their parameters, tracks found in the 
COT are extrapolated inward and matched to hits in the 
silicon detector. The SVX II has five layers of silicon 
strips at radial distances ranging from 2.5 cm to 10.6 cm 
from the beamline. Three of the five layers are double- 
sided planes with r — z strips oriented at 90° relative to 
r — (p strips, and the remaining two layers are double- 
sided planes with strips oriented at ±1.2° angles rela- 
tive to the r — (f) strips. The SVX II detector consists of 
three longitudinal barrels, each 29 cm in length, and cov- 
ers approximately 90% of the pp interaction region. The 
SVX II provides precise information on the trajectories 
of long-lived particles (decay length), which is used for 
the identification of displaced, secondary track vertices 
of B and D hadron decays. An innermost single-sided 
silicon layer (LOO), installed at 1.5 cm from the beam, 
further improves the resolution for vertex reconstruction 
j23 [. Outside of the SVX II, two additional layers of sil- 
icon assist pattern recognition and extend the sensitive 
region of the tracking detector to \ri\ « 2 [23|. These 
intermediate silicon layers (ISL) are located between the 
SVX II and the COT and consist of one layer at a radius 
of 23 cm in the central region, |?7| < 1, and two layers 
in the forward region 1 < jryj < 2, at radii of 20 and 29 
cm. The component of a charged particle's momentum 
transverse to the beam (pr) is determined with a res- 
olution of (jp^/pT ~ 0.07% pt {pt in GeV/c) for tracks 
with pt > 2 GeV/c. The excellent momentum resolution 
yields precise mass resolution for fully reconstructed B 
and D decays, which provides good signal-to-background. 
The typical resolution on the reconstructed position of 
decay vertices is approximately 30 ^m in the transverse 
direction, effective to identify vertices from charmed me- 
son decays, which are typically displaced by 250 fim from 
the beam. In the longitudinal direction, the resolution 
is approximately 70 /im, allowing suppression of back- 
grounds from charged particles originating from decays 
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of distinct heavy hadrons in the event. 



IV. ONLINE SAMPLE SELECTION 

The CDF II trigger system plays an important role 
in this measurement. Identification of hadronic decays 
of heavy-flavored mesons is challenging in the Tevatron 
collider environment due to the large inelastic pp cross 
section and high particle multiplicities at 1.96 TeV. In 
order to collect these events, the trigger system must re- 
ject more than 99.99% of the collisions while retaining 
good efficiency for signal. In this Section, we describe 
the CDF II trigger system and the algorithms used in 
collecting the samples of hadronic D decays in this anal- 
ysis. 

The CDF II trigger system has a three-level architec- 
ture: the first two levels, level 1 (LI) and level 2 (L2), are 
implemented in hardware and the third, level 3 (L3), is 
implemented in software on a cluster of computers using 
reconstruction algorithms that are similar to those used 
off line. 

Using information from the COT, at LI, the ex- 
tremely fast tracker (XFT) [24] reconstructs trajectories 
of charged particles in the r — (j) plane for each proton- 
antiproton bunch crossing. Events are selected for fur- 
ther processing when two tracks that satisfy trigger cri- 
teria on basic variables are found. The variables include 
the product of any combination of two particles' charges 
(opposite or same sign), the opening angle of the two 
tracks in the transverse plane (Ac/)), the two particles' 
transverse momenta, and their scalar sum. 

At L2 the silicon vertex trigger (SVT) [25| incorporates 
information from the SVX II detector into the trigger 
track reconstruction. The SVT identifies tracks displaced 
from the pp interaction point, such as those that arise 
from weak decays of heavy hadrons and have sufficient 
transverse momentum. Displaced tracks are those that 
have a distance of closest approach to the beamline (im- 
pact parameter do) inconsistent with having originated 
from the pp interaction point (primary vertex). The im- 
pact parameter resolution of the SVT is approximately 
50 /xm, which includes a contribution of 35 /zm from the 
width of the pp interaction region. The trigger selections 
used in this analysis require two tracks, each with im- 
pact parameter typically greater than 120 /im and smaller 
than 1 mm. In addition, the L2 trigger requires the trans- 
verse decay length {L^y) to exceed 200 /xm, where L^y is 
calculated as the projection of the vector from the pri- 
mary vertex to the two-track vertex in the transverse 
plane along the vectorial sum of the transverse momenta 
of the tracks. The trigger based on the SVT collects 
large quantities of long-lived D hadrons, rejecting most 
of the prompt background. However, through its impact- 
parameter-based selection, the SVT trigger also biases 
the observed proper decay time distribution. This has 
important consequences for the results of this analysis, 
which will be discussed in Sec. [Xj 



TABLE II. Typical selection criteria for the three versions of 
the displaced-tracks trigger used in this analysis. The crite- 
ria refer to track pairs. The pr, do, and 77 requirements are 
applied to both tracks. The 5Dpt refers to the scalar sum of 
the pt of the two tracks. The Y2pt threshold in each of the 
three vertical portions of the table identifies the high-pr (top), 
medium-pT (middle), and low-pT (bottom) trigger selections. 



Level- 1 



Level-2 



Level-3 



PT > 2.5 GeV/c PT > 2.5 GeV/c 
Ept > 6.5 GeV/c Ept > 6.5 GeV/c 
Opposite charge Opposite charge 
A<l) < 90° 2° <A(j)< 90° 

0.12 < do < 1-0 mm 



PT > 2.5 GeV/c 
Ept > 6.5 GeV/c 
Opposite charge 
2° <A(j)< 90° 
0.1 < do < 1-0 mm 



Lxy > 200 ^im Lxy > 200 /xm 

|A2;o| < 5 cm 
H < 1-2 



PT > 2 GeV/c PT > 2 GeV/c 
Ept > 5.5 GeV/c Ept > 5.5 GeV/c 
Opposite charge Opposite charge 
A(j) < 90° 2° <A(j>< 90° 

0.12 < do < 1.0 mm 
Lxy > 200 fim 



PT > 2 GeV/c 
Ept > 5.5 GeV/c 
Opposite charge 
2° <A(j)< 90° 
0.1 < do < 1.0 mm 
L^^y > 200 fim 
\Azo\ < 5 cm 
\V\ < 1-2 



PT > 2 GeV/c PT > 2 GeV/c 
EpT>4GeV/c EPT>4Ge^ 



_ GeV/c 
A<j) < 90° 2° <A(j>< 90° 

0.1 < do < 1.0 mm 
Lxy > 200 fim 



PT > 2 GeV/c 
Ept > 4 GeV/c 
2° <A(j)< 90° 
0.1 < do < 1.0 mm 
Lxy > 200 fim 
\Azo\ < 5 cm 
hi < 1-2 



The L3 trigger uses a full reconstruction of the event 
with all detector information, but uses a simpler track- 
ing algorithm and preliminary calibrations relative to the 
ones used off line. The L3 trigger retests the criteria im- 
posed by the L2 trigger. In addition, the difference in z 
of the two tracks at the point of minimum distance from 
the primary vertex, A^o, is required not to exceed 5 cm, 
removing events where the pair of tracks originate from 
different collisions within the same crossing of p and p 
bunches. 

Over the course of a single continuous period of Teva- 
tron collisions (a store), the available trigger bandwidth 
varies because trigger rates fall as instantaneous lumi- 
nosity falls. Higher trigger rates at high luminosity arise 
from both a larger rate for real physics processes as well 
as multiplicity-dependent backgrounds in multiple pp in- 
teractions. To fully exploit the available trigger band- 
width, we employ three main variants of the displaced- 
tracks trigger. The three selections are summarized in 
Table [H] and are referred to as the low-p^, medium-p^, 
and high-pT selections according to their requirements 
on minimum transverse momentum. At high luminosity, 
the higher purity but less efficient high-pr selection is 
employed. As the luminosity decreases over the course 
of a store, trigger bandwidth becomes available and the 
other selections are utilized to fill the available trigger 
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bandwidth and maximize the charm yield. The rates are 
controlled by the application of a prescale, which rejects 
a predefined fraction of events accepted by each trigger 
selection, depending on the instantaneous luminosity. 



V. SUPPRESSING DETECTOR-INDUCED 
CHARGE ASYMMETRIES 

The procedure used to cancel detector-induced asym- 
metries is briefly outlined here, while a detailed mathe- 
matical treatment is given in Appendix El 

We directly measure the observed "raw" asymmetry: 



A{D°) 



jVobs(i^°)-iVobs(i^°) 
7Vobs(i?")+iVobs(^°)' 



that is, the number of observed 13° decays into the se- 
lected final state (tt+tt" or K^K^) minus the number of 
D'^ decays, divided by the sum. 



0} 

E 

E 
>, 
w 
< 




-0.2 



1.5 2^ 
p^(7i3) [GeV/c] 



FIG. 1. Observed asymmetry between the number of recon- 
structed _D*+ and D*~ mesons as a function of the soft pion's 
transverse momentum for pure samples of D*+ -)■ L>°(-> 
7r"'"7r~)7r^ and D*~ — s- Tf{'^ n'^-K~)n^ decays. The soft 
pion transverse momentum spectrum is also shown. 

The main experimental difficulty of this measurement 
comes from the small differences in the detection effi- 
ciencies of tracks of opposite charge which may lead, if 
not properly taken into account, to spuriously-measured 
charge asymmetries. Relevant instrumental effects in- 
clude differences in interaction cross sections with matter 
between positive and negative low-momentum hadrons 
and the geometry of the main tracking system. The 
drift chamber layout is intrinsically charge asymmetric 



because of a « 35° tilt angle between the cell orienta- 
tion and the radial direction, designed to partially cor- 
rect for the Lorentz angle in the charge drift direction 
caused by crossed electric and magnetic fields. In the 
COT, different detection efficiencies are expected for pos- 
itive and negative low-momentum tracks (especially, in 
our case, for soft pious), which induce an instrumental 
asymmetry in the number of reconstructed Z3*-tagged 
and 13° mesons. Other possible asymmetries may 
originate in slightly different performance between pos- 
itive and negative tracks in pattern-reconstruction and 
track-fitting algorithms. The combined effect of these 
is a net asymmetry in the range of a few percent, as 
shown in Fig. [T] This must be corrected to better than 
one per mil to match the expected statistical precision 
of the present measurement. In order to cancel detec- 
tor effects, we extract the value of Acp{D'^ h^h^) 
using a fully data-driven method, based on an appropri- 
ate combination of charge-asymmetries observed in three 
different event samples: £'*-tagged — > h^h^ decays 
(or simply hh*), I3*-tagged D" K^Tr+ decays (Kit*), 
and untagged D" — )■ ivr~7r+ decays {Ktt). We assume 
the involved physical and instrumental asymmetries to 
be small, as indicated by previous measurements. Ne- 
glecting terms of order AcpS and 6"^, the observed asym- 
metries in the three samples are 



\hh' 



A{hh*) ^ Acp{hh) + S{TrJ' 
A{Ktt*) = AcpiKn) + (5(7r,)^'^* - 
A{Ktt) = AcpiKn) + 6{Kn)^'', 



S{Ktt 



(5) 



where 5{tts)^^ is the instrumental asymmetry for recon- 
structing a positive or negative soft pion associated with 
a h^h~ charm decay induced by charge- asymmetric in- 
teraction cross section and reconstruction efficiency for 
low transverse momentum pious; 6{tTs)^'^ is the same as 
above for tagged K^tt" and K~Tr~^ decays; and S{Kn)^'^ 
and 5{K'k)^'^ are the instrumental asymmetries for re- 
constructing a K^iT~ or a K~it^ decay for the untagged 
and the tagged case, respectively. All the above effects 
can vary as functions of a number of kinematic variables 
or environmental conditions in the detector. If the kine- 
matic distributions of soft pious are consistent in Ktt* 
and hh* samples, and if the distributions of decay 
products are consistent in Kt:* and Kt: samples, then 
SiTTsf^'' ~ (5(7r^)^^' and S{K7r)^''' « S{K7r)^''. The 
CP-violating asymmetries then become accessible as 



Acp{hh) = Aihh*) - A{Kt:*) + AiKTi). 



(6) 



This formula relies on cancellations based on two assump- 
tions. At the Tevatron, charm and anticharm mesons are 
expected to be created in almost equal numbers. Since 
the overwhelming majority of them are produced by CP- 
conserving strong interactions, and the pp initial state is 
CP symmetric, any small difference between the abun- 
dance of charm and anti-charm flavor is constrained to 
be antisymmetric in pseudorapidity. As a consequence. 
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we assume that the net effect of any possible charge asym- 
metry in the production cancels out, as long as the dis- 
tribution of the decays in the sample used for this anal- 
ysis is symmetric in pseudorapidity. An upper limit to 
any possible residual effect is evaluated as part of the 
study of systematic uncertainties (Sec. lIXp . The second 
assumption is that the detection efficiency for the D* 
can be expressed as the product of the efficiency for the 
soft pion times the efficiency for the Z?" final state. This 
assumption has been tested (Sec. IIXI) . and any residual 
effect included in the systematic uncertainties. 

Before applying this technique to data, we show that 
our approach achieves the goal of suppressing detector 
induced asymmetries down to the per mil level using the 
full Monte Carlo simulation (Appendix [BJ . The simu- 
lation contains only charmed signal decays. The effects 
of the underlying event and multiple interactions are not 
simulated. We apply the method to samples simulated 
with a wide range of physical and detector asymmetries 
to verify that the cancellation works. The simulation is 
used here only to test the validity of the technique; all 
final results are derived from data only, with no direct 
input from simulation. 



VI. ANALYSIS EVENT SELECTION 

The offline selection is designed to retain the maxi- 
mum number of — > h^h'^ decays with accurately 
measured momenta and decay vertices. Any require- 
ments that may induce asymmetries between the number 
of selected D° and mesons are avoided. The recon- 
struction is based solely on tracking, disregarding any 
information on particle identification. Candidate decays 
are reconstructed using only track pairs compatible with 
having fired the trigger. Standard quality criteria on 
the minimum number of associated silicon-detector and 
drift-chamber hits are applied to each track to ensure 
precisely measured momenta and decay vertices in three- 
dimensions [IBl- Each final-state particle is required to 
have pt > 2.2 GeV/c, \ri\ < 1, and impact parameter 
between 0.1 and 1 mm. The reconstruction of D'^ candi- 
dates considers all pairs of oppositely-charged particles in 
an event, which are arbitrarily assigned the charged pion 
mass. The two tracks are constrained to originate from 
a common vertex by a kinematic fit subject to standard 
quality requirements. The tt+tt^ mass of candidates is 
required to be in the range 1.8 to 2.4 GeV/c^, to retain 
all signals of interest and sideband regions sufficiently 
wide to study backgrounds. The two tracks are required 
to have an aziniuthal separation 2° < Ac/) < 90°, and 
correspond to a scalar sum of the two particles' trans- 
verse momenta greater than 4.5 GeV/c. We require L^y 
to exceed 200 ^m to reduce background from decays of 
hadrons that don't contain heavy quarks. We also require 
the impact parameter of the candidate with respect 
to the beam, do{D^), to be smaller than 100 fj,m to reduce 
the contribution from charmed mesons produced in long- 



lived B decays (secondary charm). In the rare (0.04%) 
occurrence that multiple 13° — ?> h^h'~ decays sharing the 
same tracks are reconstructed in the event, we retain the 
one having the best vertex fit quality. 

Figure [2] shows the K~Tr'^ mass distribution for the 
resulting sample, which is referred to as "untagged" in 
the following since no D* decay reconstruction has been 
imposed at this stage. The distribution of a sample of 
simulated inclusive charmed decays is also shown for com- 
parison. Only a single charmed meson decay per event 
is simulated without the underlying event. In both dis- 
tributions the kaon (pion) mass is arbitrarily assigned to 
the negative (positive) particle. The prominent narrow 
signal is dominated by D'^ — K^tt^ decays. A broader 
structure, also centered on the known mass, are — 
K'^iT~ candidates reconstructed with swapped K and tt 
mass assignments to the decay products. Approximately 
29 million D" and mesons decaying into K^tt^ final 
states are reconstructed. The two smaller enhancements 
at lower and higher masses than the 13° signal are due 
to mis-reconstructed — > K'^K~ and — >■ ■k^'k~ 
decays, respectively. Two sources of background con- 
tribute. A component of random track pairs that ac- 
cidentally meet the selection requirements (combinato- 
rial background) is most visible at masses higher than 2 
GeV/c^, but populates almost uniformly the whole mass 
range. A large shoulder due to mis-reconstructed multi- 
body charm decays peaks at a mass of approximately 1.6 
GeV/c2. 

In the "tagged" -samples reconstruction, we form 
D'^'^t candidates by associating with each 
candidate all tracks present in the same event. The ad- 
ditional particle is required to satisfy basic quality re- 
quirements for the numbers of associated silicon and drift 
chamber hits, to be central (I77I < 1), and to have trans- 
verse momentum greater that 400 MeV/c. We assume 
this particle to be a pion ("soft pion") and we match 
its trajectory to the D'^ vertex with simple requirements 
on relative separation: impact parameter smaller than 
600 /im and longitudinal distance from the primary ver- 
tex smaller than 1.5 cm. Since the impact parameter 
of the low-energy pion has degraded resolution with re- 
spect to those of the D'^ tracks, no real benefit is pro- 
vided by a full three-track vertex fit for the D* candi- 
date. We retain D* candidates with D'^tTs mass smaller 
than 2.02 GeV/c^. In the 2% of cases in which multiple 
D* candidates are associated with a single 13° candidate, 
we randomly choose only one D* candidate for further 
analysis. 

The D'^tTs mass is calculated using the vector sum of 
the momenta of the three particles as D* momentum, 
and the known mass in the determination of the D* 
energy. This quantity has the same resolution advantages 
of the more customary M(/i+ ft, -'^tTs) — M(ft.+/i( )~) mass 
difference, and has the additional advantage that it is in- 
dependent of the mass assigned to the decay products. 
Therefore aU D*+ D°(^ ft+/i(')")7r+ modes have the 
same D^tTs mass distribution, which is not true for the 
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Invariant K'Ti'^-mass [GeV/c^] Invariant K"7i''-mass [GeV/c^] 



FIG. 2. Comparison between the K 7r'*'-niass distributions of (a) the untagged sample and of (b) a simulated sample of 
inclusive charm decays. See text for explanation of contributions. 



mass difference distribution. 

In each tagged sample [D^ — > tt+tt^ , K^K^ 
and D''^ — > K^TT^) we require the corresponding two- 
body mass to he within 24 MeV/c^ of the known 
mass 3, as shown in Figs. [3] (a)-(c). Figures [3] (d)- 
(f) show the resulting D'^tTs mass distribution. A clean 
D* signal is visible superimposed on background com- 
ponents that are different in each channel. As will 
be shown in Sec. IVIIIl the backgrounds in the D^tTs dis- 
tributions for 13° — > TT+TT" and D'^ — >■ K'^K~ decays 
are mainly due to associations of random pions with 
real 13° candidates. In the — > K'^K~ case, there 
is also a substantial contribution from mis-reconstructed 
multi-body charged and neutral charmed decays (mainly 
Z3*+ — )• )• /C~7r+7r'')7r+ where the neutral pion is 

not reconstructed) that yield a broader enhancement un- 
derneath the signal peak. We reconstruct approximately 
215 000 D*-tagged D° tt+tt" decays, 476 000 D*~ 
tagged — >• K^K'^ decays, and 5 million £)*~tagged 
D° — > vr+A' ^ decays. 



VII. KINEMATIC DISTRIBUTIONS 
EQUALIZATION 

Because detector-induced asymmetries depend on 
kinematic properties, the asymmetry cancellation is real- 
ized accurately only if the kinematic distributions across 
the three samples are the same. Although the samples 
have been selected using the same requirements, small 
kinematic differences between decay channels may per- 
sist due to the different masses involved. We extensively 
search for any such residual effect across several kine- 



matic distributions and reweig ht the tagged D° h+h~ 
and untagged K^ir^ distributions to reproduce 

the tagged 13° — > K^tt^ distributions when necessary. 
For each channel, identical reweighting functions are used 
for charm and anti-charm decays. 

We define appropriate sideband regions according to 
the specific features of each tagged sample (Fig. [3] (a)- 
(c)). Then we compare background-subtracted distribu- 
tions for tagged /i+Zi^ )~ decays, studying a large set of 
TTs kinematic variables {pt, Vi 4'i ^.q, and zq) [26|. We 
observe small discrepancies only in the transverse mo- 
mentum and pseudorapidity distributions as shown in 
Fig. [5] (a)-(d). The ratio between the two distributions 
is used to extract a smooth curve used as a candidate- 
specific weight. A similar study of D° distributions for 
tagged and untagged D° — >• Ar~7r+ decays shows discrep- 
ancies only in the distributions of transverse momentum 
and pseudorapidity (Fig. |4|) which are reweighted accord- 
ingly. 

Background is not subtracted from the distributions 
of the untagged sample. We simply select decays with 
K~^Tr~ or ii'~7r"'"-mass within 24 MeV/c^ fr om the known 
mass, corresponding approximately to a cross-shaped 
±3(7 range in the two-dimensional distribution (Fig. [SJ. 
The background contamination in this region is about 
6%. This contamination has a small effect on the final re- 
sult. The observed asymmetries show a small dependence 
on the momentum, because detector-induced charge 
asymmetries are tiny at transverse momenta greater than 
2.2 GeV/c, as required for the D° decay products. There- 
fore any small imperfection in the reweighting of momen- 
tum spectra between tagged and untagged sample has a 
limited impact, if any. However, a systematic uncertainty 
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FIG. 3. Distributions of (a) tt+tt^, (b) K~ , and (c) Kit mass. Regions used to define the tagged samples are shaded. 
Distribution of D^\s mass for tagged (d) D'^ tv'^tt~ , (e) D" — >■ K~ , and (f) — > /(""tt^ samples selected in the shaded 
regions. 



is assessed for the possible effects of non-subtracted back- 
grounds (see Sec. lIXp . All entries in distributions shown 
in the remainder of this paper are reweighted according 
to the transverse momentum and pseudorapidity of the 
corresponding candidates unless otherwise stated. 



anti-charm, although a few parameters are determined by 
the fit independently in the two samples. The functional 
form of the mass shape for all signals is extracted from 
simulation and the values of its parameters adjusted for 
the data. The effect of this adjustment is discussed in 
Sec. IIXI where a systematic uncertainty is also assessed. 



VIII. DETERMINATION OF OBSERVED 
ASYMMETRIES 



A. Fit of tagged samples 



The asymmetries between observed numbers of Z?" and 
Z?" signal candidates are determined with fits of the D* 
(tagged samples) and (untagged sample) mass dis- 
tributions. The mass resolution of the CDF tracker is 
sufficient to separate the different decay modes of inter- 
est. Backgrounds are modeled and included in the fits. 
In all cases we use a joint binned fit that minimizes a 
combined quantity, defined as Xtot — X+ + X- j where 

and X- are the individual for the 1?° and 
distributions. Because we use copious samples, an un- 
binned likelihood fit would imply a substantially larger 
computational load without a significant improvement in 
statistical resolution. The functional form that describes 
the mass shape is assumed to be the same for charm and 



We extract the asymmetry of tagged samples by fitting 
the numbers of reconstructed D*^ events in the D^irf 

and D it J mass distribution. Because all I?" — > h^h'^ 
modes have the same D^ir'^ mass distribution, we use 
a single shape to fit all tagged signals. We also assume 
that the shape of the background from random pions as- 
sociated with a real neutral charm particle are the same. 
Systematic uncertainties due to variations in the shapes 
are discussed later in Sec. IIXI 

The general features of the signal distribution are ex- 
tracted from simulated samples. The model is adjusted 
and finalized in a fit of the D^tTs mass of copious and pure 
tagged K~Tr^ decays. We fit the average histogram of 
the charm and anti-charm samples, m = (to+ -I- m_)/2. 
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FIG. 4. Comparison between normalized kinematic distributions of the various tagged and untagged samples used in the 
analysis: (a), (c) soft pion transverse momentum, and (b),(d) pseudorapidity of hh* and Kti* events; (e) transverse 
momentum and (f) pseudorapidity of Kir and Kir* events. Tagged distributions are background-subtracted. 



where rrij^ is the D*^ mass distribution and m_. the D*" 
one. The resulting signal shape is then used in the joint 
fit to measure the asymmetry between charm and anti- 
charm signal yields. The signal is described by a Johnson 
function 27] (all functions properly normalized in the ap- 
propriate fit range), 



g-i[7 + 5 si„h-i(S^) 



J(a;|^,cr,(5,7) = 




that accounts for the asymmetric tail of the distribution, 
plus two Gaussians, J^(a;|/i, u), for the central bulk: 

psig,{m\9sig) ^JjJ{m\mD- + nj,aj,6j,'yj) + (1 - /,/) 
X [/Gl^(™|nT-D' + ^J■Gl,o■Gl) 
+ (1 - /Gl)^(™|nT-D- + MG2,CTG2)]- 

The signal parameters Ogig include the relative fractions 



between the Johnson and the Gaussian components; the 
shift from the nominal D*^ mass of the Johnson distri- 
bution's core, /ij, and the two Gaussians, Hgi{2)] the 
widths of the Johnson distribution's core, aj, and the 
two Gaussians, (7gi(2)] ^^nd the parameters Sj and 7j, 
which determine the asymmetry in the Johnson distribu- 
tion's tails. For the random pion background we use an 
empirical shape form, 

Pbkg('7i|6lbkg) = ^(to|toco +m^,6bkg,Cbkg), 

with ^(x|a, 5, c) — (x — a)''e^'^^^~''-' extracted from data 
by forming an artificial random combination made of a 
well-reconstructed D'^ meson from each event combined 
with pions from all other events. The total function used 
in this initial fit is 

iVsigpsig(m|6'sig) -I- Afbkgpbkg(nT-|6'bkg)- 
Each fit function is defined only above the threshold value 





Invariant KV-mass [GeV/c^] 



FIG. 5. Distribution of K~ n -mass as a function of K tv"- 
mass for the untagged sample. Note the logarithmic scale on 
2 axis. 



of mjjo + m^. Figure [S] shows the resulting fit which is 
used to determine the shape parameters for subsequent 
asymmetry fits. All parameters are free to float in the 
fit. 

We then fix the signal parametrization and simultane- 
ously fit the D'-'tTs mass distributions of and D*~ 
candidates with independent normalizations to extract 
the asymmetry. The parameter Sj varies independently 
for charm and anti-charm decays. The background shape 
parameters are common in the two samples and are de- 
termined by the fit. Figures [7] (a) and (b) show the 
projections of this simultaneous fit on the D^tTs mass 
distribution, for the tagged — > K~tt~^ sample. Fig- 
ures [7] (c) shows the projection on the asymmetry dis- 
tribution as a function of the D^tTs mass. The asym- 
metry distribution is constructed by evaluating bin-by- 
bin the difference and sum of the distributions in mass 
for charm ) and anti-charm (m_) decays to ob- 

tain A = (m+ — mJ)/{m-^- + to_). The variation of 
the asymmetry as a function of mass indicates whether 
backgrounds with asymmetries different from the signal 
are present. As shown by the difference plots at the bot- 
tom of Fig. [71 the fits correctly describe the asymmetry 
across the whole mass range. 

We allowed independent 5j parameters in the charm 
and anti-charm samples because the D^tTs mass distribu- 
tion for _D*+ candidates has slightly higher tails and a dif- 
ferent width than the corresponding distribution for D*^ 
candidates. The relative difference between the resulting 
5 J values does not exceed 0.5%. However, by allowing the 
parameter 5j to vary independently the x^/ndf value im- 
proves from 414/306 to 385/304. We do not expect the 



FIG. 6. Distribution of D tts mass of tagged D° K--K+ de- 
cays with fit results overlaid. The total fit projection (blue) is 
shown along with the double Gaussian bulk (dotted line) , the 
Johnson tail (dashed line) and the background (full hatching). 

source of this difference to be asymmetric background 
because the difference is maximally visible in the signal 
region, where the kinematic correlation between D'^tts 
mass and tTs transverse momentum is stronger. Indeed, 
small differences between D*~^ and D*^ shapes may be 
expected because the drift chamber has different reso- 
lutions for positive and negative low momentum parti- 
cles. Independent Sj parameters provide a significantly 
improved description of the asymmetry as a function of 
D'^tTs mass in the signal region (Fig. [7] (c)). In Sec. IIX Dl 
we report a systematic uncertainty associated with this 
assumption. No significant improvement in fit quality is 
observed when leaving other signal shape parameters free 
to vary independently for D*"*" and D*~ candidates. 

The plots in Fig. [5] show the fit results for tagged 
Z?" TT+TT" and Z?" K~K+ samples. In the 
D'^ — > K^K^ fit we include an additional component 
from mis-reconstructed multibody decays. Because sig- 
nal plus random pion shapes are fixed to those obtained 
by fitting the tagged K-k sample (Fig. [7]), the shape of 
this additional multibody component is conveniently ex- 
tracted from the combined fit to data and is described 
by 

Pmbd(f7l|0mbd) =/mbd>/("l|w_D. -|- /imbd, f mbd, ^mbd, 7mbd) 
+ (1 - /mbd)^(m-|m-£l« + "T-TT, fembd, Cmbd)- 
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FIG. 7. Results of the combined fit of tlie tagged K samples. Distribution of D'^tts mass for (a) charm, and (b) 

anti-charm decays, and (c) asymmetry as a function of the mass. Fit results are overlaid. 



The total function used to fit the KK* sample is then 

^sigpsig (to I ^sig ) + A^bkg pbkg (to I ^bkg) + A'"mbdpmbd (to | ^mbd ) • 

We observe the following asymmetries in the three 
tagged samples: 

Ai-KTT*) = (-1.86 ±0.23)%, 
A{KK*) = (-2.32 ± 0.21)%, (7) 
AiKu*) = (-2.910 ±0.049)%. 

B. Fit of the untagged sample 

In untagged Ktt decays no soft pion is associated with 
the neutral charm meson to form a D* candidate so there 
is no identification of its charm or anti-charm content. 
We infer the flavor of the neutral charm meson on a sta- 
tistical basis using the mass resolution of the tracker and 
the quasi-flavor-specific nature of neutral charm decays 
into Ktt final states. The role of mass resolution is ev- 
ident in Fig. [51 which shows the distribution of K~-k^ 
mass as a function of K^tt^ mass for the sample of un- 
tagged D° — >■ h~^h'~ decays. The cross-shaped structure 
at the center of the plot is dominated by Kn decays. 
In each mass projection the narrow component of the 
structure is due to decays where the chosen Ktt assign- 
ment is correct. The broader component is due to decays 
where the Ktt assignment is swapped. In the momen- 
tum range of interest, the observed widths of these two 
components differ by roughly an order of magnitude. Be- 
cause of the CKM hierarchy of couplings, approximately 
99.6% of neutral charm decays into a K^tt^ final state 
arc from Cabibbo- favored decays of D'^ mesons, with only 
0.4% from the doubly-suppressed decays of D mesons, 
and vice versa for A'+tt^ decays. Therefore, the narrow 



(broad) component in the AT^tt^ projection is dominated 
by D'^ (-0°) decays. Similarly, the narrow (broad) com- 
ponent in the A'+tt^ projection is dominated by D (1?°) 
decays. 

We extract the asymmetry between charm and anti- 
charm decays in the untagged sample from a simultane- 
ous binned fit of the AT+tt" and Ar^Tr"*" mass distributions 
in two independent subsamples. We randomly divide the 
untagged sample into two independent subsamples, equal 
in size, whose events were collected in the same data- 
taking period ("odd" and "even" sample). We arbitrarily 
choose to reconstruct the Ar~7r"'" mass for candidates of 
the odd sample and the AT+tt^ mass for candidates of the 
even sample. In the odd sample the decay D'^ Ar~7r+ 
is considered "right sign" (RS) because it is reconstructed 
with proper mass assignment. In the even sample it is 
considered a "wrong sign" (WS) decay, since it is re- 
constructed with swapped mass assignment. The oppo- 
site holds for the D — > AT+tt^ decay. The shapes used 
in the fit are the same for odd and even samples. The 
fit determines the number of — Ar~7r+ (RS decays) 

from the odd sample and the number of Z?" — > A'+tt^ 
(RS decays) from the even sample thus determining the 
asymmetry. Wc split the total untagged sample in half 
to avoid the need to account for correlations. The reduc- 
tion in statistical power has little practical effect since 
half of the untagged ATtt decays are still 30 (67) times 
more abundant than the tagged AT+AT^ (tt+tt") decays, 
and the corresponding statistical uncertainty gives a neg- 
ligible contribution to the uncertainty of the final result. 

The mass shapes used in the combined fit of the un- 
tagged sample are extracted from simulated events and 
adjusted by fitting the ATtt mass distribution in data. All 
functions described in the following are properly normal- 
ized when used in fits. The mass line shape of right-sign 
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FIG. 8. Results of the combined fit of the tagged — >■ 7r"'"7r and — >■ K'^ K samples. Distribution of D^tts mass for (a), 
(d) charm and (b), (e) anti-charm decays, and (c), (f) asymmetry as a function of the mass (c,f). Fit results are overlaid. 



decays is parametrized using the following analytical ex- 
pression: 

PRs("i|^Rs) =/bulk[/l^(m|m£,o +(5i,cri) 

+ (l-/i)^(m|mz,o-f ,^2,^2)] 
+ (1 - /buik)=^(TO|5, c, moo + (5i), 

where 

Sr{m\h, c, /i) = e''(™-^)Erfc(c(m - ^)), 

with Erfc(a;) — J^°° e~*^ dt. We use the sum of 

two Gaussians to parametrize the bulk of the distribu- 
tion. The function .^{ni] h, c, /i) describes the lower-mass 
tail due to the soft photon emission. The parameter /bulk 
is the relative contribution of the double Gaussian. The 
parameter /i is the fraction of dominant Gaussian, rel- 
ative to the sum of the two Gaussians. The parame- 
ters ^1(2) are possible shifts in mass from the known 



mass [3|. Because the soft photon emission makes the 
mass distribution asymmetric, the means of the Gaus- 
sians cannot be assumed to be the same. Therefore m jjo 
is fixed in the parametrization while ^1(2) are determined 
by the fit. The mass distribution of wrong-sign decays, 
Pws("^; ^'ws)j is parametrized using the same functional 
form used to model RS decays. The mass distribution 
of 13° — >■ 7r+7r~ decays is modeled using the following 
functional form: 

p^^{m\9^^) ==/buik[/i^("^|"^o + 

(l-/i)^(m|mo-H52,^2)] 
+ /ti5^(m|6i,ci,mi) 

+ (1 - /bulk - /tl)=3^(™|62,C2,™2)- 

The bulk of the distribution is described by two Gaus- 
sians. Two tail functions ^(m; b, c, /x) are added for the 
low- and high-mass tails due to soft photon emission and 
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incorrect mass assignraent, respectively. The shifts in 
mass, <^i(2)j from the empirical value of the mass of tttt 
decays assigned the Ktt mass, mg ~ 1.96736 GeV/c^, 
are free to vary. The mass distributions of the partially 
reconstructed multibody charm decays and combinato- 
rial background are modeled using decreasing exponen- 
tial functions with coefficients 5mbd and 6comb, respec- 
tively. 

The function used in the fit is then 

^RSPRS("T-|^RS) + ^WSPWS(™|^WS) 



NcomhPcomhim\k 



comb } 



where A'^rs, TVws, ^tttt, N^hd, ^comb are the event yields 
for right-sign decays, wrong-sign decays, D° — ?> tt+tt" 
decays, partially reconstructed decays, and combinatorial 
background, respectively. 
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FIG. 9. Average (m) of the distribution of K'^-k~ mass in 
the even sample and K~tt'^ mass in the odd sample with fit 
projections overlaid. 

The mass is fit in the range 1.8 < m < 2.4 GeV/c^ 
to avoid the need for modeling most of the partially re- 
constructed charm meson decays. The ratio iVR,s/.^mbd 
and the parameter 6mbd are fixed from simulated inclu- 
sive and decays. The contamination from par- 
tially reconstructed decays is negligible for masses 
greater that 1.8 GeV/c^. The result of the fit to the 
distribution averaged between odd and even samples is 
shown in Fig. O In this preliminary fit we let vary the 



number of events in each of the various components, the 
parameters of the two Gaussians describing the bulk of 
the — >■ h^h'~ distributions, and the slope of the com- 
binatorial background 6comb- We assume that the small 
tails are described accurately enough by the simulation. 
This preliminary fit is used to extract all shape parame- 
ters that will be fixed in the subsequent combined fit for 
the asymmetry. 

Odd and even samples are fitted simultaneously us- 
ing the same shapes for each component to determine 
the asymmetry of RS decays. Because no asymmetry 
in — > 7r"''7r~ decays and combinatorial background is 
expected by construction, we include the following con- 



straints: Ntnr 



and 7V+„^, 



N. 



comb ■ 



The param- 



eters TVrs, N+js, TV^g, N+^^ and iV^^d are de- 

termined by the fit independently in the even and odd 
samples. Figures [TUKa) and (b) show the fit projections 
for odd and even samples. Figure [TU] (c) shows the pro- 
jection of the simultaneous fit on the asymmetry as a 
function of the Kt: mass. The observed asymmetry for 
the D" K-n+ RS decays is 



A{KTr) 



-0.832 ±0.033)%. 



(8) 



IX. SYSTEMATIC UNCERTAINTIES 



The measurement strategy is designed to suppress sys- 
tematic uncertainties. However, we consider a few resid- 
ual sources that can impact the results: approxima- 
tions in the suppression of detector-induced asymme- 
tries; production asymmetries; contamination from sec- 
ondary D mesons; assumptions and approximations in 
fits, which include specific choice of analytic shapes, dif- 
ferences between distributions associated with charm and 
anti-charm decays, and contamination from unaccounted 
backgrounds; and, finally, assumptions and limitations of 
kinematic reweighting. 

Most of the systematic uncertainties are evaluated by 
modifying the fit functions to include systematic varia- 
tions and repeating the fits to data. The differences be- 
tween results of modified fits and the central one are used 
as systematic uncertainties. This usually overestimates 
the observed size of systematic uncertainties, which in- 
clude an additional statistical component. However, the 
additional uncertainty is negligible, given the size of the 
event samples involved. Sources of systematic uncer- 
tainty are detailed below. A summary of the most sig- 
nificant uncertainties is given in Table IIIII 



A. Approximations in the suppression of 
detector-induced effects 

We check the reliability of the cancellation of all 
detector-induced asymmetries on simulated samples as 
described in Appendix |B] The analysis is repeated on 
several statistical ensembles in which we introduce known 
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FIG. 10. Results of the combined fit of the untagged — >■ K -k^ sample. Distribution of D^tts mass for (a) charm, and (b) 
anti-charm decays, and (c) asymmetry as a function of the mass. Fit results are overlaid. 



CP-violating asymmetries in the D'^ — > h^h!^ decays 
and instrumental effects (asymmetric reconstruction effi- 
ciency for positive and negative soft pions and kaons) de- 
pendent on a number of kinematic variables (e.g., trans- 
verse momentum). These studies constrain the size of 
residual instrumental effects that might not be fully can- 
celled by our method of linear subtraction of asymme- 
tries. They also assess the impact of possible correlations 
between reconstruction efficiencies of D° decay-products 
and the soft pion, which are assumed negligible in the 
analysis. We further check this assumption on data by 
searching for any variation of the observed asymmetry as 
a function of the proximity between the soft pion and the 
charm meson trajectories. No variation is found. 

Using the results obtained with realistic values for the 
simulated effects, we assess a AAcp{hh) = 0.009% un- 
certainty. This corresponds to the maximum shift, in- 
creased by one standard deviation, observed in the re- 
sults, for true CP- violating asymmetries in input ranging 
from -5% to +5%. 



B. Production asymmetries 

Charm production in high-energy pp collisions is dom- 
inated by CP-conserving cc production through the 
strong interaction. No production asymmetries are ex- 
pected by integrating over the whole phase space. How- 
ever, the CDF acceptance covers a limited region of the 
phase space, where CP conservation may not be ex- 
actly realized. Correlations with the pp initial state 
may induce pseudorapidity-dependent asymmetries be- 
tween the number of produced charm and anti-charm (or 
positive- and negative-charged) mesons. These asymme- 
tries are constrained by CP conservation to change sign 
for opposite values of 77. The net effect is expected to 



vanish if the pseudorapidity distribution of the sample is 
symmetric. 

To set an upper limit to the possible effect of small 
residual rj asymmetries of the samples used in this anal- 
ysis, we repeat the fits enforcing a perfect 77 symmetry 
by reweighting. We observe variations of lS.Acp{KK) = 
0.03% and A^cpItttt) = 0.04% between the fit resuhs 
obtained with and without re-weighting. We take these 
small differences as an estimate of the size of possible 
residual effects. The cancellation of production asymme- 
tries achieved in pp collisions (an initial CP-symmetric 
state) recorded with a polar-symmetric detector provide 
a significant advantage in high-precision CP-violation 
measurements over experiments conducted in pp colli- 
sions. 



C. Contamination of D mesons from B decays 

A contamination of charm mesons produced in h- 
hadron decays could bias the results. Violation of CP 
symmetry in 6-hadron decays may result in asymmet- 
ric production of charm and anti-charm mesons. This 
may be large for a single exclusive mode, but the effect 
is expected to vanish for inclusive B D^X decays [28j. 
However, we use the impact parameter distribution of Z)" 
mesons to statistically separate primary and secondary 
mesons and assign a systematic uncertainty. Here, by 
"secondary" we mean any originating from the decay 
of any b hadron regardless of the particular decay chain 
involved. In particular we do not distinguish whether the 
meson is coming from a D*^ or not. 

If /b is the fraction of secondary Z?*^ mesons in a given 
sample, the corresponding observed asymmetry A can 
be written as a linear combination of the asymmetries 
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for primary and secondary mesons: 

A = fBA{D° secondary) + (1 - fB)A{D" primary). (9) 

The asymmetry observed for secondary D'^ mesons can be 
expressed, to first order, as the sum of the asymmetry you 
would observe for a primary D'^ sample, plus a possible 
CP-violating asymmetry in inclusive B — >■ D'^X decays, 

A{D" sec.) = AcpiB D°X) + A{D° prim.). (10) 

Hence, combining Eq. Q and Eq. (|10p. the asymmetry 
observed in each sample is given by 

A = fsAcpiB D"X) + A{D" primary). (11) 

Because the fraction of secondary Z?" mesons is inde- 
pendent of their decay mode, we assume fsi'^'^*) — 
fB{KK*) = fBiKn*). The contribution of CP violation 
in 6-hadron decays to the final asymmetries is written as 

A{hh) = fB{K'!i)Acp{B D"X) + AcpiD° hh), 

(12) 

where /b is estimated in the untagged K 7r+ sample be- 
cause the two terms arising from the tagged components 
cancel in the subtraction provided by Eq. ([5]). In this 
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FIG. 11. Impact parameter distribution of D° candidates in 
the — >■ K~n^ signal region. Top plot with data and fit 
projections overlaid uses a logarithmic scale vertically. Bot- 
tom plot shows fractional difference between data and the fit 
on a linear scale. 

analysis, the contamination from secondary decays is 



reduced by requiring the impact parameter of the can- 
didate, do{D^), not to exceed 100 ^m. The fraction Jb of 
residual mesons originating from B decays has been 
determined by fitting the distribution of the impact pa- 
rameter of untagged D'^ — >■ K~tt'^ decays selected within 
±24 MeV/c^ of the known D° mass Q. We use two 
Gaussian distributions to model the narrow peak from 
primary D'^ mesons and a binned histogram, extracted 
from a simulated sample of inclusive B D'^X decays, 
to model the secondary component. Figure [TT] shows the 
data with the fit projection overlaid. A residual con- 
tamination of 16.6% of -B — >■ D^X decays with impact 
parameter lower than 100 fira is estimated. To constrain 
the size of the effect from Acp{B D'^X) we repeat 
the analysis inverting the impact parameter selection, 
namely requiring do{D'^) > 100 /xm. This selects an al- 
most pure sample of D'^ K^tt'^ decays from B decays 
[Jb — 1). We reconstruct about 900 000 decays with an 
asymmetry, A(if7r) = (— 0.647±0.172)%, consistent with 
(—0.832 ± 0.033)%, the value used in our measurement. 
Using Eq. (ITUl) we write the difference between the above 
asymmetry and the asymmetry observed in the central 
analysis (Eq. ^^), A{dQ > 100 ^m) - A{do < 100 /im), 
as 

(1 - fB)Acp{B D°X) = (-0.18 ± 0.17)%. (13) 

Using Jb = 16.6% we obtain Acp(,B D^X) = 
(—0.21 ± 0.20)% showing that no evidence for a bias 
induced by secondary mesons is present. Based on 
Eq. , we assign a conservative systematic uncertainty 
evaluated as /b^cp(S ^ DX) = /b/(1 - /b)A = 
0.034%, where /s equals 16.6% and A corresponds to the 
0.17% standard deviation of the difference in Eq. (IT^ . 



D. Assumptions in the fits of tagged samples 

1. Shapes of fit functions 

The mass shape extracted from simulation has been 
adjusted using data for a more accurate description of 
the observed signal shape. A systematic uncertainty is 
associated with the finite accuracy of this tuning and 
covers the effect of possible mis-modeling of the shapes 
of the fit components. 

Figure [12] shows a comparison between the shape ex- 
tracted from the simulation and the templates used in 
the fit after the tuning. It also shows an additional tem- 
plate, named "anti-tuned" , where the corrections that 
adjust the simulation to data have been inverted. If 
/(m) is the template tuned on data, and g{'m) is the 
template extracted from the simulation, the anti-tuned 
template is constructed as h{m) = 2/(m) — g{m). We 
repeat the measurement using the templates extracted 
from the simulation without any tuning, and those cor- 
responding to the anti-tuning. The maximum variations 
from the central fit results, A^c'p(7r+7r") = 0.009% 
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FIG. 12. Shape of D tts mass as extracted from simulation 
without tuning, with data tuning and with anti-data tuning. 

and AAcp{K~^K~) — 0.058%, are assigned as system- 
atic uncertainties. The larger effect observed in the 
_D" — > K^K~ case comes from the additional degrees 
of freedom introduced in the fit by the multibody-decays 
component. 

In addition, we perform a cross-check of the shape used 
for the background of real mesons associated with 
random tracks. In the analysis, the shape parameters of 
— > h^h^ fits are constrained to the values obtained in 
the higher-statistics tagged — ?> K~tt~^ sample. If the 
parameters are left fioating in the fit, only a negligible 
variation on the final result (< 0.003%) is observed. 



2. Charge-dependent mass distributions 

We observe small differences between distributions of 
D'^TTs mass for positive and negative D* candidates. 
These are ascribed to possible differences in tracking res- 
olutions between low-momentum positive and negative 
particles. Such differences may impact our results at 
first order and would not be corrected by our subtrac- 
tion method. To determine a systematic uncertainty, 
we repeat the fit in several configurations where various 
combinations of signal and background parameters are 
independently determined for positive and negative D* 
candidates. The largest effects are observed by leaving 
the background shapes to vary independently and con- 
straining t he p arameter 6j of the Johnson function to be 
the same [26|. The values of the shape parameters in 
— h^h~ fits are always fixed to the ones obtained 
from the — )■ if "tt"*" sample. The maximum variations 
with respect to the central fits, A^c'p(7r+7r^) — 0.088% 



and AAcp{K~^K ) = 0.027%, are used as systematic 
uncertainties. 



3. Asymmetnes from residual backgrounds 

A further source of systematic uncertainty is the ap- 
proximations used in the subtraction of physics back- 
grounds. In the K^K^ sample we fit any residual back- 
ground contribution, hence this uncertainty is absorbed 
in the statistical one. However, in the tt+tt" and K~Tr~^ 
cases we assume the residual backgrounds to be negli- 
gible. Using simulation we estimate that a 0.22% and 
0.77% contamination from physics backgrounds enters 
the ±24 MeV/c^ tt+tt" and K~tt^ signal range, respec- 
tively. The contamination in the tt+tt^ sample is domi- 
nated by the high mass tail of the Z?" K^n^ signal. 
The asymmetry of this contamination is determined from 
a fit of the tagged K~tt^ sample. The contamination of 
the K^t:^ sample is dominated by the tail from partially 
reconstructed D'~' decays. The fit of the tagged K~^K~ 
sample provides an estimate of the asymmetry of this 
contamination. In both cases we assign a systematic un- 
certainty that is the product of the contaminating frac- 
tion times the additional asymmetry of the contaminant. 
This yields a maximum effect of 0.005% on the measured 
asymmetries for both Z?" tt+tt^ and D'^ — >■ K~^K~ 
cases. 



E. Assumptions in the fits of untagged samples 

1. Shapes of fit functions 

We follow the same strategy used for the tagged case 
to assign the systematic uncertainty associated with pos- 
sible mis-modeling of the shapes in fits of the untagged 
sample. 

Figure [T3] shows the comparison between templates ex- 
tracted from the simulation without any tuning, those 
tuned to data (and used in the central fit) , and the anti- 
tuned ones. We repeat the fit using the templates from 
simulation and the anti-tuned ones. The maximum vari- 
ation from the central fit, AA{Ktt) — 0.005%, is used as 
the systematic uncertainty. 

2. Charge-dependent mass distributions 

In the untagged case we expect the mass shapes of all 
components to be the same for charm and anti-charm 
samples. However, we repeat the simultaneous fit un- 
der different assumptions to assign the systematic uncer- 
tainty associated with possible residual differences. The 
parameters of the Gaussian distributions used to model 
the bulk of the mass distributions are left free to vary in- 
dependently for the charm and anti-charm samples, and 
separately for the right-sign, wrong-sign, and D — ^ tt+tt" 
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FIG. 13. Shapes of if^Tr^ mass from simulation without tuning, with data tuning, and with anti-data tuning for (a) right-sign 
and (b) wrong-sign K^n^ decays, and for (c) tt+tt" decays. 



components. We assume no difference between mass dis- 
tributions of combinatorial background and partially re- 
constructed decays. The differences between estimated 
shape parameters in charm and anti-charm samples do 
not exceed 3cr, showing compatibility between the shapes. 
A systematic uncertainty of 0.044% is obtained by sum- 
ming in quadrature the shifts from the central values of 
the estimated asymmetries in the three different cases. 



3. Asymmetries from residual physics backgrounds 



In the measurement of the asymmetry of Cabibbo- 
favored — ?► K^tt^ decays, we neglect the contribution 
from the small, but irreducible, component of doubly- 
Cabibbo-suppressed (DCS) 1?° — > K^tt^ decays. Large 
CP violation in DCS decays may bias the charge asymme- 
try we attribute to 13° — )■ X~7r+ decays. We assign a sys- 
tematic uncertainty corresponding to fDCS^CpiD^ 
K+TT-) = /dcsA = 0.013%, where foes = 0.39% is 
the known Q fraction of DCS decays with respect to 
Cabibbo- favored decays and A — 2.2% corresponds to 
one standard deviation of the current measured limit on 
the CP-violating asymmetry Acp{D'^ — > K'^tt') as re- 
ported in Rcf. [3|. 

In the central fit for the untagged 13" K^n^ sample, 
no asymmetry in D'^ tt^tt" decays or combinatorial 
background is included, as expected by the way the un- 
tagged sample is defined. We confirm the validity of this 
choice by fitting the asymmetry with independent param- 
eters for these two shapes in the charm and anti-charm 
samples. The result corresponds to a AA{KTr) — 0.011% 
variation from the central fit. 



F. Limitations of kinematic reweighting 



The tagged event samples are reweighted after sub- 
tracting the background, sampled in signal mass side- 
bands. We constrain the size of possible residual system- 
atic uncertainties by repeating the fit of tagged — > 
h^h^ after a reweighting without any sideband subtrac- 
tion. The variation in observed asymmetries is found to 
be negligible with respect to other systematic uncertain- 
ties. 

In reweighting the untagged sample we do not sub- 
tract the background. The signal distributions are ex- 
tracted by selecting a mass region corresponding ap- 
proximately to a cross-shaped window of ±3(7 in the 
two-dimensional space {M{K'^Tr~), M{K^Tr'^)). To as- 
sign a systematic uncertainty we extract the signal dis- 
tributions and reweight the data using a smaller cross- 
shaped region of ±2a (i.e. within 16 MeV/c^ from the 
nominal D'^ mass). The background contamination de- 
creases from 6% to 4%. We repeat the analysis and find 
A{KTr) = (— 0.831±0.033)% corresponding to a variation 
from the central fit of < 0.001%, thus negligible with re- 
spect to other systematic uncertainties. 



G. Total systematic uncertainty 



Tabic mil summarizes the most significant systematic 
uncertainties considered in the measurement. Assuming 
them independent and summing in quadrature, we ob- 
tain a total systematic uncertainty of 0.11% on the ob- 
served CP-violating asymmetry of D'^ — > tt+tt^ decays 
and 0.09% in D° — >■ K~^K~ decays. Their sizes are ap- 
proximately half of the statistical uncertainties. 
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TABLE III. Summary of most significant systematic uncertainties. The uncertainties reported for the last three sources result 
from the sum in quadrature of the contributions in the tagged and untagged fits. 



Source 


Acp{tt+-k ) 


:%] Acp{K+K-) [%] 


Approximations in the suppression of detector-induced effects 


U.UUa 


u.uuy 


Production asymmetries 


0.040 


0.030 


Contamination of secondary D mesons 


0.034 


0.034 


Shapes assumed in fits 


0.010 


0.058 


Charge-dependent mass distributions 


0.098 


0.052 


Asymmetries from residual backgrounds 


0.014 


0.014 


Limitations of sample reweighting 


< 0.001 


< 0.001 


Total 


0.113 


0.092 



X. FINAL RESULT 

Using the observed asymmetries from Eqs. ([7]) and ([5]) 
in the relationships of Eq. ([5]), we determine the time- 
integrated CP-violating asymmetries in tt+tt" and 
D" K+R- decays to be 

Acp(7r+7r") = (+0.22 ±0.24 (stat) ± 0.11 (syst))% 
AcpiK+R-) = (-0.24 ±0.22 (stat) ± 0.09 (syst))%, 

corresponding to CP conservation in the time-evolution 
of these decays. These are the most precise determina- 
tions of these quantities to date, and significantly im- 
prove the world's average values. The results are also in 
agreement with theory predictions pol - ls^ . 

A useful comparison with results from other experi- 
ments is achieved by expressing the observed asymmetry 
as a linear combination (Eq. (j?])) of a direct component, 
App, and an indirect component, A™p, through a coeffi- 
cient that is the mean proper decay time of charm mesons 
in the data sample. The direct component corresponds 
to a difference in width between charm and anti-charm 
decays into the same final state. The indirect component 
is due to the probability for a charm meson to oscillate 
into an anti-charm meson being different from the prob- 
ability for an anti-charm meson to oscillate into a charm 
meson. 

The decay time of each D*^ meson, t, is determined as 

_ Lxy _ TO_DO 

c(P7)t cpT 

where {P^)t = Pt/tt^d<^ is the transverse Lorentz factor. 
This is an unbiased estimate of the actual decay time 
only for primary charmed mesons. For secondary charm, 
the decay time of the parent B meson should be sub- 
tracted. The mean decay times of our signals are deter- 
mined from a fit to the proper decay time distribution of 
sideband-subtracted tagged decays fFig. [T4t . The fit in- 
cludes components for primary and secondary D mesons, 
whose shapes are modeled from simulation. The simula- 
tion is used to extract the information on the mean de- 
cay time of secondary charmed decays, using the known 
true decay time. The proportions between primary and 
secondary are also determined from this fit and are con- 
sistent with results of the fit to the impact parameter 



in data (Sec. IIX Cp . We determine a mean decay time 
of 2.40 ± 0.03 and 2.65 ± 0.03, in units of D'^ hfetime, 
for — >■ 7r"'"7r~ and K^K^ decays, respectively. 

The uncertainty is the sum in quadrature of statistical 
and systematic contributions. The small difference in the 
two samples is caused by the slightly different kinematic 
distributions of the two decays, which impacts their trig- 
ger acceptance. 

Each of our measurements defines a band in the 
{A'^% A^^'p) plane with slope - {t) /r (Eq. (g])). The same 
holds for BAiAR and Belle measurements, with slope 
— 1 [13, [Hi, due to unbiased acceptance in decay time. 
The results of this measurement and the most recent B- 
factories' results are shown in Fig. [T5l which displays 
their relationship. The bands represent ±1(t uncertain- 
ties and show that all measurements are compatible with 
CP conservation (origin in the two-dimensional plane). 
The results of the three experiments can be combined as- 
suming Gaussian uncertainties. We construct combined 
confidence regions in the (^™p, ^cp) pl^ne, denoted with 
68% and 95% confidence level ellipses. The correspond- 
ing values for the asymmetries are A^p(I?° — )■ tt+tt^) — 
(0.04 ± 0.69)%, ^'cp(-D° ^ TT+TT") = (0.08 ± 0.34)%, 
Afp{D^ ^ K+K-) = (-0.24±0.41)%, and A'^${D° ^ 
K+K^) = (0.00±0.20)%, in which the uncertainties rep- 
resent one-dimensional 68% confidence level intervals. 



A. CP violation from mixing only 

Assuming negligible direct CP violation in both decay 
modes, the observed asymmetry is only due to mixing, 
Acp{h+h-) « A'Sf> {t}/T, yielding 

A'gl>{n+TT~) ^ (+0.09 ±0.10 (stat) ±0.05 (syst))% 
A'l^'^{K+K^) = (-0.09 ±0.08 (stat) ± 0.03 (syst))%. 

Assuming that no large weak phases from non-SM contri- 
butions appear in the decay amplitudes, is indepen- 
dent of the final state. Therefore the two measurements 
can be averaged, assuming correlated systematic uncer- 
tainties, to obtain a precise determination of CP violation 
in charm mixing: 

A'gl>{D^) = (-0.01 ± 0.06 (stat) ± 0.04 (syst))%. 



21 



15 



10 



xlO'' 



(a) 



D**^ D" (^7l+7t ) 7t+ + c.c. 

X^/ndf = 87/59 

• Data (5.9 fb-'; 

Fit 

I Non-prompt component 




0"^ 


1 

2 


4 


6 8 

t/x 



















CM 

O 

CD 
Q. 

</} 
CD 
CS 
"P 
T3 

CO 

O 



30 



20 



10 



xlO'' 



(b) 




D** ^ D" K*K ) 7t+ + c.c. 

X^/ndf = 59/59 



• Data (5.9 fb"') 

Fit 

I Non-prompt component 




FIG. 14. Distribution of proper decay time (in units of Z?" lifetime) for sideband-subtracted tagged (a) — !> tt+tt and (b) 
— > K~ data. Fit results are overlaid including the component from secondary charmed mesons (red). 



This corresponds to the following upper limits on CP 
violation in charm mixing: 

\A'^${D^)\ < 0.13 (0.16)% at the 90 (95)% C.L. 

The bias toward longer-lived decays of the CDF sam- 
ple offers a significant advantage over i?-factories in sen- 
sitivity to the time-dependent component, as shown in 
Figs. [111(a), (c). 

B. Direct CP violation only 

Assuming that CP symmetry is conserved in charm 
mixing, our results are readily comparable to mea- 
surements obtained at _B-factories; Acp{t^^t^^) = 
(0.43 ± 0.52 (stat) ± 0.12 (syst))% and Acp{K+K~) = 
(-0.43 ± 0.30 (stat) ± 0.11 (syst))% from Belle, and 
Acp{TT+Tr-) = (-0.24 ± 0.52 (stat) ± 0.22 (syst))% and 
Acp{K+K-) = (0.00 ± 0.34 (stat) ± 0.13 (syst))% from 
BMAR (Figs.[in](b)-(d)). The CDF result is the world's 
most precise. 

C. Difference of asymmetries 

A useful comparison with theory predictions is 
achieved by calculating the difference between the asym- 
metries observed in the _D" — > K~^K~ and — > tt+tt" 



decays (AAcp)- Since the difference in decay-time accep- 
tance is small, A{t)/T — 0.26 ±0.01, most of the indirect 
CF-violating asymmetry cancels in the subtraction, as- 
suming that no large CF-violating phases from non-SM 
contributions enter the decay amplitudes. Hence AAcp 
approximates the difference in direct CP- violating asym- 
metries of the two decays. Using the observed asymme- 
tries from Eq. ([7]), we determine 

AAcp ^Acp{K+K-) - Acp{tt+^-) 
^AAtp + A^^%A{t)lT 
=A{KK*)- A{t:t:*) 
= (-0.46 ±0.31 (stat) ±0.12 (syst))%. 

The systematic uncertainty is dominated by the 0.12% 
uncertainty from the shapes assumed in the mass fits, 
and their possible dependence on the charge of the D* 
meson. This is determined by combining the differ- 
ence of shifts observed in Sees. IIX D II and IIX D 21 in- 
cluding correlations: (0.058 - 0.009)% = 0.049% and 
(-0.027- 0.088)% = 0.115%. Smaller contributions in- 
clude a 0.009% from the finite precision associated to the 
suppression of detector-induced effects (Sec. IIX Ap . and 
a 0.005% due to the 0.22% background we ignore under 
the 1?° TT+TT" signal (Sec lIXD3l) . The effects of pro- 
duction asymmetries and contamination from secondary 
charm decays cancel in the difference. 
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FIG. 15. Comparison of the present results with Belle and BA3AR measurements of time-integrated CP-violating asymmetry 
in (a) — >■ tt^tt" and (b) D" — K'^ K~ decays displayed in the {A^qp, A'qp) plane. The point with error bars denotes the 
central value of the combination of the three measurements with one-dimensional 68% confidence level uncertainties. 



We see no evidence of a difference in CP violation be- 
tween K^K'^ and — > tt+tt^ decays. Figure 
[T7] shows the difference in direct asymmetry (Ayl^p) as 
a function of the indirect asymmetry compared with ex- 
perimental results from BjBAR and Belle O, [IH- The 
bands represent ilcr uncertainties. The measurements, 
combined assuming Gaussian uncertainties, provide 68% 
and 95% confidence level regions in the (AA^p, A^^p) 
plane, denoted with ellipses. The corresponding val- 
ues for the asymmetries are Aj4^p = (—0.37 ± 0.45)%, 
(-0.35 ±2.15)%. 



XI. SUMMARY 



the charmed mesons used in these measurements are 
2.40 ± 0.03 units of D° hfetime in the D° tt+tt" 
sample and 2.65 ± 0.03 units of lifetime in the 
K^K^ sample. Assuming negligible CP violation 
in — 7r+7r~ and — K~^K~ decay widths (direct 
CP violation) , the above results, combined with the high- 
valued average proper decay time of the charmed mesons 
in our sample, provide a stringent general constraint 
on CP violation in D° mixing, \A'gf>{D°)\ < 0.13% at 
the 90% confidence level. The results probe significant 
regions of the parameter space of charm phenomenology 
where discrimination between SM and non-SM dynamics 
becomes possible [s^ [s!] . 



In summary, we report the results of the most sensitive 
search for CP violation in singly-Cabibbo-suppressed 
I?" — >■ 7r+7r~ and — )■ K~^K~ decays. We reconstruct 
signals of 0(10^) -D*-tagged decays in an event sam- 
ple of pp collision data corresponding to approximately 
5.9 fb~^ of integrated luminosity collected by a trigger 
on displaced tracks. A fully data-driven method to can- 
cel instrumental effects provides effective suppression of 
systematic uncertainties to the 0.1% level, approximately 
half the magnitude of the statistical uncertainties. 

We find no evidence of CP viola- 
tion and measure Acp{D'^ — > tt+tt^) = 
(+0.22 ± 0.24 (stat) ± 0.11 (syst))% and 
Acp{D° K+K-) = (-0.24 ± 0.22 (stat) ± 0.09 (syst))% 
These are the most precise determinations from a single 
experiment to date^ and supersede the correspond- 
ing results of Ref. [131. The average decay times of 
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FIG. 16. Comparison of the present results with results from Belle and BA3AR assuming (a), (c) no direct, or (b), (d) no 
indirect CP violation. In each plot the la band of the B-factories' average is displayed in blue, while the new average that 
includes the CDF result is shown in green. 
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Appendix A: Method to suppress detector 
asymmetries 

A mathematical derivation of the concepts described 
in Sec. |V] follows. We measure the CP- violating asym- 
metry by determining the asymmetry between number 



of detected particles of opposite charm content A — 
{N+-N^)/_{N++N^), where N+ and N_ are the number 
of Z?" and decays found in three different data sam- 
ples: £'*-tagged h^h^ decays (or simply hh*), 
£'*-tagged — K~Tr^ decays (Kir*) and untagged 
D° — > Ar~7r+ decays (Kit). We show that the combi- 
nation of asymmetries measured in these three samples 
yields an unbiased estimate of the physical value of Acp 
with a high degree of suppression of systematic uncertain- 
ties coming from detector asymmetries. In the discussion 
we always refer to the true values of kinematic variables 
of particles. The measured quantities, affected by ex- 
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A^^ [%] 

FIG. 17. Difference between direct CP-violating asymmetries 
in the K'^ K~ and tt^tt" final states as a function of the indi- 
rect asymmetry. Belle and B^AR measurements are also re- 
ported for comparison. The point with error bars denotes the 
central value of the combination of the three measurements 
with one-dimensional 68% confidence level uncertainties. 

perimental uncertainties, play no role here since we are 
only interested in counting particles and all detection ef- 
ficiencies are assumed to be dependent on true quantities 
only. 

1. 73*-tagged D° h+h~ 



and h~^h~; is the branching fraction of 

£)*+ D°iT+ and D*~ — > D°it~, assumed to be charge- 
symmetric; Ehh is the detection efficiency of the h^h~ 
pair from the 1?° decay; and 2^+ and Eg- are the detec- 
tion efficiencies of the positive and negative soft pion, 
respectively. Conservation of four-momenta is implic- 
itly assumed in all densities. Densities are normalized as 
/dp*p*±(p*) = 1 = J dpsdph+dph-phh'{Ph+,Ph-,Ps\p*) 
for each . The difference between event yields is there- 
fore 

f 

N+ - iV- = — Bh^ dp^dpsdp,^+dp^- 

X Phh*{Ph+,Ph-,Ps \p*)£hh{Ph+,Ph-) 

X {p^+{p^)B+^es+iPs) - p*^{p*)B^,^es-{ps)} 

iV, f 

= — B*D^ dp^ dpsdph+ dph -ehhiPh+,Ph-) 

X Phh*iPh+,Ph-iPs \p*)p*{p*)Bhhes{ps) 
X [(1 + Sp4p,)) (1 + Acp) (1 + SesiPs)) 
- (1 - Sp,{p,)) (1 - Acp) (1 - SssiPs))], 

where we have defined the following additional quantities: 
P* = (1/2) (p*+ + P*-), Sp* = - + p*-), 

S,;, - (1/2) (S+ + Acp ^ Acp{hh) - (5+ - 

BHh)/iBth + B^h)^ = (l/2)(e.+ + Es-), and fc, = 
(es+—Ss-){es++es-)- Expanding the products we obtain 

iV+ - 7V_ ^N^B*^^Bhh j dp^dpsdpi^+dph- p*{p*)e s{ps) 

X Phh-{Ph+,Ph-,Ps \p*)£hh{Ph+,Ph-) 
X [Acp + Sp^{p^) + 6es{ps) 
+ Acp6p,:{p,:)6esips)]- 



Assuming factorization of efficiencies for reconstruct- 
ing the neutral charmed meson and the soft pion, we 
write 

N± =—B*jj^ dp^dpsdph+dph- p*±{p*)Bi^^ 

X Phh'{Ph+,Ph-,Ps \ P*)ehh{Ph+ ,Ph-)£s±{Ps), 

where N* is the total number of D*^ and D*^ mesons; 
Pi,,Ps,Ph+ ,Ph- are the three-momenta of the D* , soft tt, 

and h~ , respectively; and p*_ are the densi- 
ties in phase space of D*~^ and D*~ mesons (function 
of the production cross sections and experimental ac- 
ceptances and efhciencies) ; phh* is the density in phase 
space of the soft pion and h^h" pair from 13° decay; 

and 5^,^ are the branching fractions of h^h^ 



Since the CP symmetry of the pp initial state ensures 
that (5p*(ph,) = — p*), the second and fourth term 
in brackets vanish when integrated over a domain sym- 
metric in r]. In a similar way we obtain 

N+ + =N^B*jj^Bhh Jdp^dpsdph+dph- p*{p*)SsiPs) 

X Phh'{Ph+,Ph-,Ps \p*)£hh{Ph+,Ph-) 
X [1 + AcpSes{ps)AcpSp*{p*) 
+ Ses{ps)Sp*ip*)]. 

The second term in brackets is small with respect to Acp 
and can be neglected, while the third and fourth terms 
vanish once integrated over a domain symmetric in rj. 
Hence the observed asymmetry is written as 



25 



A{hh*) = 



= Acp{h+h )+ / dp.Nl^ (p^)fes(Ps), where 



N+-N_\ 
N+ + N_) 

Jdp^dph+ dph- p* {p*)phh- {Ph+ ,Ph-,Ps\ P*)£hh {ph+ ,Ph-)£s jPs) 
Jdp^ dph+ dph- dpsp* {p*)phh' {Ph+ ,Ph-,Ps\p* )£hh {Ph+ ,Ph-)ss{Ps) 

I 



(Al) 
(A2) 



is the normahzed density in phase space of the soft pion 
for the events included in our sample. 



2. D*-tagged D° K'n'' 



Assuming factorization of efficiencies for reconstruct- 
ing the neutral charmed meson and the soft pion, we 
write 



N± =-^B*^^ j dp^Apsdp^dpK P*±{p*)B 

X PKn''{PK,Pn,Ps \ P*)£KTn±{PK ,Ptv)£s±{Ps), 



where p,r and pk are the three- momenta of the pion and 

kaon, p*j^^ is the density in phase space of the soft pion 
and Kn pair from the decay, and B]^^ are the 
branching fractions of — )• K~'k~^ and — >■ /C+tt", 
and SK--IT+ and Sk+h- arc the detection efficiencies of 
the K~Tr'^ and K'^Tr~ pairs from D'^ and D° decay. The 
difference between charm and anti-charm event yields is 



written as 

-/V* f 

N+ - iV_ =-^B*jj^ dp^dpsdp^dpKPK7r*{PK,Pn,Ps \p*) 

X [p^+{p^)B+^eK-'K+{PK,P'K)Ss+iPs) 
- p^_{p^)B-^eK+n-iPK,Pn)Ss-iPs)] 

=-^B*jj^Bkt, / dp^dpsdp^dpKP*{p*)es{Ps) 

X PK7r-iPK,P7v,Ps \P*)£KniPK,P-K) 

X {{1 + Sp4p,)){1 + Acp) 

X (1 + feif7r(PK,P7r))(l + SSsiPs)) 

- {l-dp,{p,))il-Acp) 

X (1 - fex7r(PK,P7r))(l - fes(Ps))}, 

where we have defined the following additional quantities: 

Bk. = {l/2){B+_^ + B^J, AcP ^ Acp{K-k) = {B+^ - 
BK-K)/iB^^ + Bj^^), skt: = (l/2)(eK-7^+ + £K+^-), and 
5£ktt = (£/f-7r+ - £K+T7-)/{eK-Tv+ + e/f+Tr-)- Expand- 
ing the products and observing that all terms in Sp^{p^) 
vanish upon integration over a symmetric p* domain, we 
obtain 

N+ - =Nt^B*j^^BK-K jdp^dpsdp^dpK p*ip*)£s{Ps) 

X PKt,'{PK,P7t,Ps \p*)£ktt{pk,P,t) 

X {AcP + SeK7r{PK,P7v) + S£s{Ps) + •••}, 

where we have neglected one term of order AcpS^. Sim- 
ilarly, 

--N^BI)^Bk7t Jdp^dpsdp,,dpKP*{p*)£s{Ps) 

X Pk-k'{pk,P7,,Ps \p*)£k-k{pk,P7,) 
X [1 + Acp5eK,,{pK,PT,) + AcpSssiPs) 

+ 5SK7r{PK,P7r)S£s{Ps)]- 

If we neglect all terms of order AcpS and 6"^, we finally 
obtain 



N+ + N_ 



A{Ktt 



N+ - N_ 



Ktt* 



N+ + N. 
where h^l'{pK,Pn) 



=Acp{K-TT+)+ J dp^h^l' {pK,Prr)S£KniPK,P7:) + J dp, /if (p, (p, ) , (A3) 

_ Jdp^dp,p^(p^)pK7T'{PK,P-K,Ps \p*)£k-k{pk,P-k)£s{Ps) 



J dp*dpyrdpKdpsP*{p*)pK7T'iPK,P7r,Ps \p*)£K7:{PK,P7T)£siPs)' 



(A4) 
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and h^'^* {ps) (the Kir analogous to h^^' {p^) in Ea. lX2|) 
are the normahzed densities in phase space of 'k,K and 
soft TT, respectively, for the events included in our sample. 



3. Untagged K'-r^ 



In this case 



=^ / dpQdpj,dpK Po±{po)B^^ 



N+ - =^Bk-k j dpodp^dpK 

X Po{Po)PK'niPK,P7v \Po)eK7r{PK,PTi) 
X {(1 + d>o(po))(l + Acp){l + SeKnipK,P^)) 
- (1 - 5po{po)){l - Acp){l - SeK7r{pK,PTT))} 

where we have defined the following quantities po = 
(1/2) {po+ + PQ-) and Spo = (po+ - Po-)/(po+ + Po-)- 
Assuming rj symmetry of the po integration region, 

N+~N-=NoBk-k J dpodpTrdpKPQ{pQ)pK7riPK,P7T \ Po) 
X £K7T{PK,P7r)[AcP - SeK^{pK,P-,T)]- 

Similarly we obtain 

N+ + iV_ ^NqBkt, J dpodp^dpK Po{po)PK7TiPK ,P7r I Po) 
X eKn{PK,PiT)[l + AcpSeK7T{PK,Pii)], 

and neglecting the second term in brackets, 



A{KTr) 



N+ - iV_ 



AcpiK-TT^ 



hKZ{PK,Pn) 



J dpoPo{po)p%^{pK,PTv I Po)eK7r{PK,Pn) 

J dpodp^dpKPoipo)p%^ipK,P^ \pa)£K-KipK,P7v) 

I 



dp^ dpK hjil [pK ,P7r)5eK7r [PK , Ptt ) , where 



is the normalized density in phase space of the Kt: system 
in the events included in our sample. 



4. Combining the asymmetries 

By combining the asymmetries measured in the three 
event samples we obtain 

I 



(A5) 



A{hh*)-A{Kn*) + A{K^) = Acp{h+h-) + J dp./if * (p,)fe,(p,) 

- Acp{K^TT+) - J dpKdpyrh^l {pK,P7T)SeK-^{PK,P7r) - J dpshf"" {ps)5es{ps) 

+ Acp{K-n+) + J dpKdp.,h'^l{pK,Pn)SeK^{pK,P^) = Acp{h+h'), (A6) 



where we assumed h^'^ {ps) = h!l^ {ps), and 
h§Z' (pkjPtv) = h^ZiPK,Piv)- The last two equalities 
are enforced by appropriate kinematic reweighing of the 
event samples. We need to equalize distributions with re- 



spect to the true momenta while we only access the distri- 
butions with respect to the measured momenta. Hence 
the assumption that event samples that have the same 
distribution with respect to the measured quantities also 
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[GeV/c] 

FIG. 18. Curves corresponding to simulated ratios of effi- 
ciencies for reconstructing positive versus negative pions as a 
function of transverse momentum. 



have the same distribution with respect to the true quan- 
tities is needed. 

The mathematical derivation shows that for small 
enough physics and detector-induced asymmetries, the 
linear combination of the observed asymmetries used in 
this measurement achieves an accurate cancellation of the 
instrumental effects with minimal impact on systematic 
uncertainties. 



Appendix B: Monte Carlo test of the analysis 
technique 

We tested the suppression of instrumental effects by re- 
peating the analysis in simulated samples in which known 
instrumental and physics asymmetries were introduced. 
Many different configurations for the input asymmetries 
were tested, covering a rather extended range, to ensure 
the reliability of the method independently of their ac- 
tual size in our data. For each configuration, 0(10^) 
decays were simulated to reach the desired 0.1% sensitiv- 
ity. Only the 13° — tt+tt" sample was tested although 
the results are valid for the Z?" — )• K^K~ case as well. 

We test cancellation of instrumental effects arising 
from different reconstruction efficiencies between positive 
and negative particles, which in general depend on the 
particle species and momentum. Furthermore, the relia- 
bility of the suppression should not depend on the actual 
size of CP violation in — > K^tt^ and — > tt+tt^ 
decays. 

We repeated the measurement on statistical ensembles 



where the above effects are known and arbitrarily var- 
ied using a combination of event-specific weights applied 
to the true values of simulated quantities. Each ensem- 
ble consists of approximately one thousand trials. We 
compare the resulting observed asymmetry ^'^p(7r7r) to 
the one given in input, A^p°(7r"''7r~), by inspecting the 
distribution of the residual, AAcpiTTTr) = A'^p(7r+7r") — 

We first investigate the individual impact of each ef- 
fect. We scan the value of a single input parameter across 
a range that covers larger variations than expected in 
data and assume all other effects are zero. First a px- 
dependent function that represents the dependence ob- 
served in data (see Fig. [T|) is used to parametrize the 
soft pion reconstruction efficiency ratio as e(7r"'")/e(7r~) = 
Erf {1.5 ■ pt + A), where pr is in GeV/c and various val- 
ues of the constant A have been tested so that the effi- 
ciency ratio at 0.4 GeV/c spans the 0.6-1 GeV/c range 
as shown in Fig. 1181 Then, the kaon reconstruction effi- 
ciency ratio e{K^) /e{K^) is varied similarly in the 0.6-1 
GeV/c range. Finally, a range -10% < Acp < 10% 
is tested for the physical CP-violating asymmetry in 
D° — > K^TT^ and — > tt+tt^ decays. 

The results are shown in Fig. [TO] (empty dots). The 
cancellation of instrumental asymmetries is realized at 
the sub-per mil level even with input effects of size much 
larger than expected in data. 

Figure [in] (filled dots) shows the results of a more com- 
plete test in which other effects are simulated, in addi- 
tion to the quantities varied in the single input parame- 
ter scan: a pr-dependent relative efficiency e(7r+)/e(7r~), 
corresponding to 0.8 at 0.4 GeV/c, e{K-)/€{K+) = 98%, 
Acp{Ktt) = 0.8% and Acp{t^t^) = 1.1%. Larger varia- 
tions of the residual are observed with respect to the pre- 
vious case. This is expected because mixed higher-order 
terms corresponding to the product of different effects 
are not canceled and become relevant. 

Finally we tested one case with more realistic val- 
ues for the input effects. The pT dependence of 
e(7r+)/e(7r^) is extracted from fitting data (Fig.[T]) to be 
distributed as Erf (2.49 pr), with px in GeV/c. We used 
e{K+)/e{K-) w e{K+iT-)/(i{K--K+) = 1.0166, in which 
the approximation holds assuming equal efficiency for re- 
constructing positive and negative pions at pt > 2 GeV/c 
(37| . We assume AcpiKn) = 0.1%, ten times larger 
than the current experimental sensitivity. A —5% < 
Acp(7r7r) < 5% range is tested in steps of 0.5% for the 
physical asymmetry to be measured. The results are 
shown in Fig. 1201 The maximum observed bias is of 
the order of 0.02%, one order of magnitude smaller than 
the statistical resolution on the present measurement. 
The observed bias is (0.0077 ± 0.0008)% averaged over 
the A(7p(7r7r) range probed. These results, which extend 
to the K^K^ case, demonstrate the reliability of our 
method in extracting a precise and unbiased measure- 
ment of CP violation in meson decays into K~^K~ 
and 7r+7r~ final states, even in the presence of sizable 
instrumental asymmetries. 
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FIG. 19. Asymmetry residual as a lunction of the input quantity varied. Other effects are assumed zero (empty dots) or 
different from zero (filled dots). 



The results discussed in this appendix are used in nal results due to neglecting higher order terms in Eq. 
Sec. IIXI to estimate a systematic uncertainty on the fi- (jH]), including possible non-factorization of h'^h ~ and tt^ 

reconstruction efficiencies. 
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